Andreev bound states for superconducting-ferromagnetic box 
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Within the microscopic Bogoliubov-de Gennes (BdG) formalism an exact quantization condition 
for Andreev bound states of the ferromagnetic-superconducting hybrid systems of box geometry is 
derived and a semi-classical formula for the density of states is obtained. The semi-classical formula 
is shown to agree with the exact result, even when the exchange field h, is much larger than the 
superconductor order parameter, provided h is small compared with the Fermi energy. 
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Mesoscopic hybrid systems formed from ferromagnets 
(F) in contact with superconductors (S) exhibit inter- 
esting transport properties resulting from the suppres- 
sion of the electron-hole correlation in the ferromag- 
ne ^ 3 i i 2 j 2 j 4££ j Z£^ These transport phenomena are inti- 
mately related to the influence of an exchange field on 
the density of states (DOS) of clean ferromagnetic films 
in contact with superconductor, which has been investi- 
gated both experimentally^ and theoretically^. In FS 
systems, Andreev bound states below the bulk supercon- 
ducting gap are spin split by the exchange field of the 
ferromagnet. In a quasi-classical treatment of the dif- 
fusive regime it was shown that sub-gap features in the 
DOS of FS hybrids can be understood from the behavior 
of the length distribution of the classical trajectories ex- 
isting in the ferromagnetic region, which depends purely 
on the geometry and the boundaries of the ferromagnetic 
region. 

In this work we calculate the Andreev levels of an FS 
box consisting of a clean ferromagnetic region with a su- 
perconductor attached to one side, as shown in Fig. ^ 
This geometrical arrangement is a generalization of that 
investigated in Ref.Uot where the size of the system along 
the FS interface was infinite. The discrete energy spec- 
trum of the FS box system is obtained by solving the 
microscopic Bogoliubov-de Gennes (BdG) equationiii. 
We derive an exact quantization condition without using 
the frequently applied Andreev's approximation. Our ex- 
act quantum description of the FS hybrid is an extension 
of the commonly used model developed by Blonder et 
al>i£ and by Saint-James and de Gennesi 3 - for normal- 
superconducting systems. The mismatch in the effective 
masses and Fermi energies of the ferromagnet and the su- 
perconductor are taken into account in our calculations 
and the tunnel barrier at the interface is modeled using 
a Dirac delta potential. The treatment of the problem 
is based on an adaptation of a method developed in our 
previous workia to the case of FS systems. 

Starting from this exact quantization condition, we 
give a semi-classical expression for the sub-gap density 
of states (DOS) for exchange fields which are much less 
than the Fermi energy. The semi-classical DOS is ex- 
pressed in terms of the classical return probability of 
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FIG. 1: A rectangular shape of ferromagnetic dot (F) in con- 
tact with a superconductor (S). 

the electrons which depends only on the geometry of 
the F region and its boundaries. Besides the DOS, 
from our derivation an explicit expression for this return 
probability is also found. Based on the quasi-classical 
Green's function approach a similar expression has been 
derived for the DOS of ferromagnetic layers in contact 
with clean superconductor in Ref. |ul Our DOS expres- 
sion for FS systems can be regarded as an extension 
of the Bohr-Sommerfeld formula developed for normal- 
superconducting hybrids 15 . In this work we compare the 
DOS obtained from the exact quantum calculations with 
that found from our semi-classical formula. We show 
that a good agreement between the two calculations can 
be expected only for small enough exchange fields (much 
less than the Fermi energy). 

The BdG equation for the FS systems can be written 

as 
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where H = p 2 /2m e g + V(r) — fi is the single-particle 
Hamiltonian, /i = Ep F \ are the Fermi energies, 
m g = mF,TOs are the effective masses in the F/S re- 
gions, is a two-component wave function, E is the 
quasi-particle energy measured from the Fermi energy 
Ep . Here a — ±1 for spin up/down states and +/ — 
refer to the electron/hole like quasi-particle excitation. 
The tunnel barrier V(r) at the FS interface and the 
exchange energy h(r) are modeled in a usual way by 
V(x,y) = Uq5(x) and h(x,y) — h®(x), where O is the 
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unit step function. We also adopt the usual step-function 
modeliiiS for the pair potential and take A(r) = A6(i). 
It is easy to see that the Hamiltonian is separable and 
the ansatz for the wave functions in the F region can be 
written as 
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while in the S region the wave functions have the form 
^ m ,a(x 7 y) 
where 

/ \ 2 

(2c) 




l(F) 



\ 



1± 



E 



(F) 



q (±) 

Xm(y) 

7± 



k 



(S) 



\ 



l±i 



VE 2 - A 2 



(S) 



y/2/W sm(mny/W), 

±i arccos(£'/A) 



,(2d) 

(2e) 
(2f) 



Here m is a fixed integer and the Fermi wave numbers in 
the F/S regions are given by fc F = 



2m F E F F) /h 2 and 
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, j, J 2msE^ /h 2 . The wave functions satisfy the 
Dirichlet boundary conditions at the boundary of the 
F region except for the FS interface where the match- 
ing conditions 14 should be applied. The four coefficients 
a±,c± in Eqs. I|2a(l and l|2bjl are determined from these 
matching conditions. One can find the following secular 
equation for the eigenvalues E of the FS system for fixed 
mode index m and spin state a: 

Im { 7+ DW(E, h) DH(E, h)} = 0, (3a) 

where Im{.} stands for the imaginary part and 
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and Z = 2mpUo/h 2 is the normalized barrier strength. 
The number of propagating modes for the electron/hole 



are the maximum of m for which km}a is a real num- 

where M = 
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k^W/ir and [•] stands for the integer part. For non- 
propagating modes k^ a have to be replaced with their 
imaginary part and the functions sin and cos with the 
functions sinh and cosh, respectively. For fixed m and a 
the solutions of Eq. i(3*a|) for E give the discrete sub-gap 
energy spectrum (E < A). These levels are numerically 
exact, i.e., no Andreev's approximation is used. 

We now calculate density of states below the gap. In 
what follows, we assume that there is no mismatch and 
tunnel barrier at the FS interface (mp = ms , E F — E F 
and Z = 0). For simplicity, we shall omit the superscript 
F and S in the wave numbers and the Fermi energies. 
In Andreev's approximation, i.e., for \E + <rh\ <C Ep 

we have k m . a w q m<a and D { „t}r(E,h) w k£,l e~ lk ™-<* d . 



Therefore, the quantization condition (|3all can be simpli- 
fied 

(k { ml - k ( m }a] d - arccos (E/A) 
I m A E ) = - =n. (4) 

The density of states for energies below the gap (\E\ < 
A ) is 
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where i? mn .o- the solutions of Eq. Q. Using J5J) the 
DOS becomes 
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where M = M^ZJ +1 , the number of propagating modes 
for spin-up-hole. Applying the Poisson summation 
formula 1 ^ to the summation over m and keeping only 
the non-oscillating term one finds 
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where the -E-dependent m* satisfies Eq. |@J for a given 

n and a. To simplify I m ,a we Taylor expand km,cr in 
terms of E + ah in first order (which is consistent with 
the Andreev's approximation) and find 
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where vp is the Fermi velocity. In our approximation, 
Mq ps M. From Q and using (JSJ) we obtain 
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and £o = hvp/A is the coherence length in the bulk su- 
perconductor. 



N(E) 




FIG. 2: Counting function N(E) as a functions of E/A ob- 
tained from the exact (solid line) and the semi-classical cal- 
culation (dashed line) for h/A = 0.1 (h/E F = 0.0025). The 
other parameters are M = 217.7, d/W = 0.7, A/E F = 0.025. 
The insets show the enlarged part of the main frame. 



be shown that P(s) is the classical probability that an 
electron entering the billiard at the FS interface returns 
to the interface after a path length s. The distribution 
P(s) is normalized to one, i.e., r o °° P{s)ds = 1. Note 
that the result for the DOS, given in Eq. (|10() differs from 
that obtained in Ref. by a factor multiplying P(s) in 
the summation. 
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FIG. 3: Counting function N(E) as a functions of E/A ob- 
tained from the exact (solid line) and the semi-classical cal- 
culation (dashed line) for h/A = 10.0 (h/E F = 0.25). The 
other parameters are the same as in Fig. [5] The insets show 
the enlarged part of the main frame. 

Using © and performing the derivatives in Q we find 
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is a purely geometry-dependent function, £ c (E) = 
Co/V 1 - E 2 /A 2 and s n , a (E) is given by Eq. ©. It can 



FIG. 4: The density of states g(E) as a functions of E/A 
(in units of 2qn) obtained from the exact (solid line) and the 
semi-classical calculation (dashed line) for parameters given 
in Fig. H 
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FIG. 5: The density of states g(E) as a functions of E/A 
(in units of 2qn) obtained from the exact (solid line) and the 
semi-classical calculation (dashed line) for parameters given 
in Fig. 01 



To compare the exact DOS obtained from the quanti- 
zation condition (|3afl with that calculated from the semi- 
classical expression <|1(J[) . we introduce the integrated 
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FIG. 6: Counting function N(E) as a functions of E/A ob- 
tained from the exact (solid line) and the semi-classical calcu- 
lation (dashed line) for h/A = 40.0 (h/E F = 1.0). The other 
parameters are the same as in Fig. 



DOS: N(E) = J dE'g(E'). From fT® we have 
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and F(s) = f° P{s')ds' = 9 (s - 2d) ^1 - 4&'7s 2 is the 
integrated length distribution. Note that F(oo) = 1 since 
P(s) is normalized to 1. 

The small parameter used in reaching the semi- 



classical result is h/Ep. Figures [5] and [3] show a com- 
parison between the semi-classical result Eq. (jl 1|> and 
the exact result Eq. ©, when h/E F = 0.0025 and 0.25, 
respectively. For better resolution, in the insets enlarged 
parts of the main frames are shown. Figures 0] and [S] 
show the corresponding densities of states (in units of 
2gN, where £>n = rn-pAj (ttH 2 ) is the DOS of the ferromag- 
netic region including spin up/down states and A = Wd 
is the area of this region) with singularities given by the 
semi-classical formula s„ !0 -(F s i ng ) = 2c?, where a = ±1 
and n is such that F s i ng < A. These figures show that 
the semi-classical formula given by Eqs. (|10|l - l|ll|l yield 
good agreement with the exact result, even for large val- 
ues of h/A, provided h/Ep is small. Figure shows a 
comparison with the exact result when the latter condi- 
tion is violated. Clearly the agreement is poor in this 
limit. 

In summary we have shown that a semi-classical treat- 
ment of the clean limit yields an expression for the DOS 
in terms of the classical return probability P(s), which 
in turn is known analytically. This formula is analogous 
to the quasi-classical result of Ref. llOl where P(s) is not 
known analytically and must be determined via a nu- 
merical simulation. We have also shown that the semi- 
classical formula agrees very well with the exact result for 
small exchange field compared with the Fermi energy. 
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